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Various orbital-dependent exchange-only potentials are studied which exhibit correct long-range asymptotic behaviour. We 
present the first application of these potentials for polymers and by one of these potentials for molecules. Kohn-Sham type 
calculations have been carried out for polyethylene in order to make valuable comparison of these potentials with each other as 
well as with Hartree-Fock and exchange-only LDA (Xa) methods. The difference between total energies and highest occupied 
orbital energies obtained with Hartree-Fock methods and with localized exchange-potentials is larger for this polymer than for 
atoms or molecules. Various properties of the band structure are also calculated. The band gap strongly depends on the basis 
set. The larger basis set makes the Kohn-Sham eigenvalue gap too low at about 4.4 eV while the minimal basis set results in 
value close to the experimental gap (~ 8.8 eV). For the total energy and the exchange energy, the various orbital-dependent 
exchange-only and Hartree-Fock results differ only slightly, but for the highest occupied orbital energy the difference is more 
pronounced. The Kohn-Sham band gap obtained with the optimized effective potetial method is corrected with the exchange 
contribution to the derivative discontinuity of the exchange-correlation potential. The corrected band gap obtained with the 
Slater's exchange potential is 9.7 eV close to the experiment. 



1. Introduction 



Density Functional Theory (DFT) is well established 
in theoretical physics and quantum chemistry 0|. The 
application of the theory to investigate the properties 
of materials has been made in the framework of the so- 
called second generation density functionals in which the 
kinetic energy is expressed in terms of orbitals but the 
exchange and correlation energy E^c in terms of the elec- 
tron density. This leads to the one-electron Kohn-Sham 
(KS) equations. In many cases E^c is treated in the lo- 
cal density approximation (LDA). However, for atoms 
and molecules gradient correction to the electron density, 
such as the generalized gradient approximation (GGA) is 
necessary in order to improve the results. 

Density functionals of the third generation treat both 
the kinetic and exchange energy exactly in terms of 
single-particle orbitals and only the correlation energy 
needs to be approximated, either in terms of the den- 
sity or in terms of the orbitals. The corresponding one- 
particle equations are known as the optimized effective 
potential method (OEP) derived for the first time quite 
some time ago by Sharp and Horton and later by Talman 
and Shadwick . These OEP integro-differential equa- 
tions yield orbitals which are implicit functionals of the 
density p since the orbitals come from a local potential. 
Therefore, one can use the following notation for such a 
class of potentials: r)}]. The Hohenberg-Kohn 

theorem Q], applied to noninteracting systems, ensures 
that the ground-state determinant, hence, all the occu- 
pied orbitals are unique functionals of the density. 

However, the OEP method is still rather complicated. 
Krieger, Li and lafrate (KLI) transformed and approxi- 
mated these equations into a manageable form and ap- 
plied them to atoms |^,^ where they yield results nearly 
identical to OEP [||. Later the KLI approach was also 
used for molecules [Q and for certain semiconductors 
but not for polymer chains. Other calculations with 
a Kohn-Sham potential which is believed to be close to 
the exact KS potential resulted in much narrower gap 
than the experiment |]l0| . 

The KLI approximation to OEP provides a self- 
interaction free exchange potential and proper — i 
asymptotic behaviour as r — > cx) Other approximate 



potentials ^^[{wi}] are reported which exhibit proper 
asymptotics as well The Koopmans theorem Q is 
satisfied by the highest occupied orbital in OEP and KLI, 
while it is violated by other density functionals such as 
LDA and GGA. For instance, while the popular GGA 
exchange-correlation energy functional [pl| provides large 
improvements in the relative energies of various systems, 
the exchange-correlation potentials and the energy eigen- 
values of the highest occupied state, are both suffering 
from significant error Also others found insignifi- 
cant improvement to semiconductor energy gap by GGA 
|l2| . In the last decades quite a number of publications 
appeared on polymer electronic structure using mainly 
Hartree-Fock ab initio methods [|l3| . Results are reported 
only the most recently, however, using DFT methods like 
LDA|,0. 

In this paper we will present results obtained by vari- 
ous orbital dependent local exchange potentials for poly- 
mers. We give the valuable comparison of different Vx 
potentials which allows one to draw definite conclusions 
about the future applicability of OEP based methods 
for periodic systems. The work presented here describes 
the calculated band gaps, band widths, band structures 
and total energy of polyethylene within the Kohn-Sham 
density functional scheme using, however, various or- 
bital dependent exchange potentials as the Slater's po- 
tential and the reasonably simplified version of the KLI 
exchange-only potential. For comparison we also give re- 
sults obtained by Hartree-Fock (HE) and simple periodic 
exchange-only LDA (Xa). 

2. Basic formalism 

We start from the following Kohn-Sham effective 
single-particle equation of ordinary density functional 
theory (for sake of simplicity we dropped the spin index 
and atomic units are used throughout the paper): 

h^V^ + Vs{r)]u^{r) = e,w,(r) (1) 

where {ui, i = 1,2, TV} are eigenfunctions of the single 
particle effective KS potential Vs{r) |]lj. i is a collective 
index for all single particle quantum numbers. 

Vs{r) ^ Vext{r) + VHir) + Vxc{r) (2) 
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As usual, Vext{j^) denotes the external and ^//(r) the 
Hartree potential. Vxc(j) is a local exchange-correlation 
potential which is formally defined as the functional 
derivative of the exchange-correlation energy 



Vxc{r, [p]) = 



Sp{v) ■ 



(3) 



Due to the Hohenberg-Kohn theorem, there exists a one- 
to-one mapping between the single particle KS potentials 
Vs{r) and densities p(r) which guarantees that the func- 
tional derivative given by Eq. (|) is defined 
The exact Exc[p] can be written as 



Exc[p] = 2 



dr'dr 



' ^^ p{r)p{r')[gH[p];r,r')-l] ^ 



(4) 

where g'^{[p]; r,r') is the pair correlation function of the 
ficticious system with interaction strength parameter < 
A < 1 and a ground state density which is independent 
of A. When employing Eq. (S) on Exc[p], 



Vxc{r,[p]) 







r - r 
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The first term of the potential is the so-called poten- 
tial part of the exchange-correlation energy mE 18 1 which 



The exact exchange part of Exc is known in terms of 
the single particle orbitals, i.e. 



is identical with Slater's potential in the exchange-only 
case, i.e. for g^ approximated by g^^'^ and which is ex- 
actly known [ pT| . To get the second, response-like term 
on the right-hand side of Eq. (5) one faces the tech- 
nical difficulty that - even for the exchange-only case - 
the functional derivative can not be determined directly. 
One way to get Uj:c(r) for a given, orbital-dependent ap- 
proximation to Exc is provided by the optimized effective 
potential method (OEP) given by Talman and Shadwick 
The starting point of the OEP method is the total 
energy functional 

EZr[P\=f:jdvu:[r) (^-iv^)^,(r) 

+ j dvp{v)vcxt{v) + \ I drdr'^^^^+Eg^^iiu,}], (6) 

In contrast to ordinary DPT, the exchange-correlation 
energy is an explicit functional of orbitals and therefore 
only an implicit functional of the density via Eq. (||) . 
The local single-particle potential appearing in Eqs. ffl) 
and (j^) is obtained by minimizing i.e. 



5Eor[{u,}] 
5vs{r) 



= 0. 



(7) 



As first pointed out by Perdew and co-workers Eq. (Q) is 
equivalent to the Hohenberg-Kohn variational principle 

mm 



w occ „ 

i?— = J drdV 



(8) 



We will not repeat the details of the OEP calculations 
here and refer the interested reader to a recent review 
of the method Also from now on we will restrict 

ourselves to the exchange-only case, i.e. we neglect cor- 
relation completely and use Eq. (|^) as approximation to 
the exchange-correlation functional. The resulting OEP 
integral equation is rather complicated to solve in prac- 
tice . In the work of Krieger and co-workers the OEP 
integral equation is analyzed and a simple approxima- 
tion is made which reduces the complexity of the original 
OEP equation significantly and at the same time keeps 
many of the essential properties of OEP unchanged 
Krieger, Li and lafrate gave an exact expression by trans- 
forming the OEP integral equations into a manageable 
form. They got the following, still exact, expression for 
vx{r) |]: 



P(r) 



-OEP —HF 



v^i) (9) 



l^ V[p.(r)Vu,(r)] 
2^ p{r) 

where (r) is the Slater potential given as the first term 
in Eq.(5) when g^ = g^=^ and pi(r) — |wi(r)p. Slater's 
potential can be written in terms of the first-order density 
matrix 7(r,r ) as follows 



^Slater, . ^ ( bi^-^)? 

" ^' 2 J p(r)|r-r' 



dr'. 



(10) 



In Eq. (9) the summation runs over the orbital index for 
all the occupied orbitals up to the highest occupied mth 
orbital (Fermi level) . The function pi is defined by 

p.(r) = ^ / dr'K°/^(r')-z;.(r')]G,(r,r')^,(r'), 
Uj(r) J 

(11) 



where Gi(r, r') is the Green's function 



ej - 



(12) 



In practical applications the last term in Eq.(9) turned 
out to be quite small in atomic systems and has small 
effect only at the atomic shell boundaries ||]. Krieger, Li 
and lafrate (KLI) || have proposed a simple approxima- 
tion where this last term is neglected completely. This 
might appear a rather crude approximation but it can 
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be interpreted as a mean-field approximation since the 
neglected terms averaged over the ground-state density 
p(r) vanish. This is not only true for finite systems 
with exponentially decaying densities but also for 
infinite systems such as solids or polymers The 
KLI-approximation to the exchange potential after some 
algebra reads 



p(r) ^ 



(13) 



-HF 



where m is the highest occupied orbital level. 



(14) 



P(r) 



dr,i,j = 1, ...,m~ 1. 



(15) 



f f„ and u^j^ are given as follows: 



v^j = / Pj(r)vf (r)dr. 



«f(r) 



^ occ 

4i: 



,<(r)u*(r>,(r')Mj(r) 



2p( 



w»(r) 



^E<(r)u,(r) J dv 



,u.(r')u*(r') 



(17) 



(18 



1 5EO^''[{u,}] 
<(r) <5u,(r) ^ 



(19) 



The term corresponding to the highest occupied orbital 
Um has been excluded from the sum in Eq. ( |l3[ ) , because 

-KLI ^ -HF rn 

^ xm xm o 

In this article we restrict ourself to different exchange- 
only methods, however, the extension of these methods 
to exchange-correlation case is in principle straightfor- 
ward . The most trivial way of accounting for Coulomb 
correlation effects is the use of a correlation energy func- 
tional either in its local or gradient corrected form. How- 
ever, recent calculations using a gradient- and orbital- 
dependent functional in combination with exact exchange 
provided excellent results for atoms only The poor 
performance of this functional for molecules must be at- 
tributed to the improper long-range component in the 
corresponding correlation hole which is needed to can- 
cel the long-range component in the exact exchange hole 



(the combined xc hole is typically short-ranged). Most 
of the approximate correlation functionals, which are de- 
rived from the homogeneous or the inhomogeneous elec- 
tron gas model, are suffering from incorrect long-range 
assymptotics and exhibit improper local behaviour [[l5[ . 
Therefore the extension of OEP based methods to in- 
clude electron correlation will probably be the subject of 
further studies in the next future 

In this paper we do not attempt to solve Eq. ( p^ ) 
for polymers but instead choose a different approach 
proposed by Gritsenko at. al ||l^. For the constants 
Wj = v^j^^ — v^f in Eq. (|l^) an alternative expression 
is proposed |]l8| in terms of orbital energies Cj . It follows 
from gauge invariance requirements, proper scaling and 
short range behaviour of the response part of Eq.(5) and 
Eq. (^3|) that w., does only depend on energy differences. 
Therefore Eq. (|lj) turns into the following much simpler 
formula: 



,,SSP 



(r) 



"x\^)+ (j.) 



(20) 



where is the Fermi level (highest occupied energy 
level). We use the notation SSP (Slater's potential -I- 
step potential) for this exchange potential. The con- 
stant 1^ is obtained from the homogeneous electron gas 
(HEG) model so that, Eq. (|2^) is exact in the HEG hmit 
(16) [|l8[. This constant is chosen as universal parameter for 
all the calculations. It is one of the main advantages 
of this expression that one can avoid the tedious ma- 
trix inversion of Eq. (|3|) and also suggests the future 
applicability of a more general class of v^c which can 
be designated by Vxci[{ui}, {ei}];r) |2^. Note that the 
summation runs over all the occupied orbitals except the 
highest one, as in Eq. ([l3|). With the step potential- 
like second term v^^^ provides a good approximation 
to the OEP exchange-potential (Eq.(9)) [|l8[ that pos- 
sesses proper short-range behaviour and the characteris- 
tic atomic-shell stepped structure. 



3. The implementation of the orbital depen- 
dent exchange potentials 

In this section we briefly give the summary of the the- 
ory for a periodic one-dimensional polymer chain based 
on gaussian lobe functions |^ . 

In the following we want to describe quasi one dimen- 
sional helical chain polymers, i.e. we consider a single 
polymer chain (extending to infinity) which may or may 
not have a helix structure (for general background on the 
LCAO method on helical polymers see, e.g. |14|). Math- 
ematically this symmetry can be expressed by a screw 
operation S in terms of a translation a along the z-axis 
combined with a rotation of angle Q about the same axis, 
i.e. formally. 
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xcosQ — ysinO 
S{a,Q)r — I xsinQ + ycosS 



(21) 



The symmetry group generated by the screw operation is 
an Abehan group and therefore its irreducible representa- 
tions are one-dimensional. Therefore the single-electron 
wavefunctions corresponding to this problem will trans- 
form according to 



Ui{k, S"r) — exp(—ikn)ui{k, r) 



(22) 



where the label k (which may be viewed as a normalized 
wavenumber) is restricted to the range —n < k < t: and 
i is simply the band index (to avoid confusion of indexes 
we note that i in the exponent is simply the imaginary 
number) . 

The electrons are assumed to doubly occupy a set of 
one-electron orbitals, Ui{k,r), of Bloch-type. These or- 
bitals are written as linear combinations of mf, atomic 
(real) basis functions. {X/j;M ~ l,---,m-fc} denotes the 
basis atomic functions in the reference unit cell and 
X]](r) ~ x^(S"r) are the corresponding basis functions 
for unit cell n. The ansatz for the orbitals then reads 



+00 

E 



nib 



c.M(^)e''"x:(r) 



(23) 



which transforms according to Eq. (p^). In this article 
the sets of index {n,n ,n ), and (/i, cr, t) are used, 
referring to cells, orbitals and contracted atomic lobe 
functions. The unknown coefficients Ci^(fc) have to be de- 
termined through the self-consistency procedure. Since 
we are dealing with an infinite chain polymer, the sum 
over the neighbouring cells extends from —00 to +00. In 
practice, however, we only take a finite number of neigh- 
bouring cells into account. In section 5. we will address 
the question of how many cells need to be taken into 
account in order to get converged ground-state energies. 

The total energy for the polymer system is expressed 
as H 

mt, -\-oc 



■Hit, -|-oc- ^ 

111/ n = — rv-i ^ 



+ {xl\e.c{[{u.}];r)\x: 



fii/ n— — oo 

+ 00 



It E 



nib +OC 



n — — 00 



n' — — 00 ^ ^ rs ' ' ^ 

lb +00 

E E p"'" ( E E p:?:^" ix 

n— — 00 

2E(4 



|r - R" 



where Zr and Rf denote the nuclear charges and coor- 
dinates within a single unit cell, R" denotes the nuclear 
coordinates in unit cell n (R" = S^Ilr). In our par- 
ticular case the orbital dependent exchange-correlation 



energy density is exci[UH}];r) = ex{[{ui}];r) and also 
Ex = J e^{[{u^}];r)dr 

According to Eq. (|) E^^^ can be given in terms of 
one-particle density matrixes as well. 



OEP 



4 ^ 



-t-00 

E 



pO,n \ ^ -pn ,■ 



X°(r)x;;(rOxg(r)x?"(rO ^ 
|r — r'l 

where xJJ, x",x^ ^'^'^ Xt the basis functions in 

which the orbitals are expanded. The upper indexes 
and n,n' ,n denote the localization of the basis function 
in the reference cell and in the infinite system, respec- 
tively. The density matrix elements for one-dimensional 
periodic systems are given by, 



occ 

E 



^c*„{k)cjrik)exp[ik{n' - n )]. (26) 



The density matrix P^i/ is computed at each iteration 
by numerical integration over the occupied part of the 
first Brillouin zone of the polymer. The interactions of 
the "reference cell" with neighbours (finite neighbouring 
interaction) are taken into account by the summation 
over the cell index n. The electron density can than be 
written as 



'°W= E EPM"x°(r)xI^(r). 



(27) 



n— — 00 /.Li^ 



The contracted gaussian-lobe Xni^) basis functions in 
the reference unit cell (22| is equal to 



1 



(28) 



where the summation runs over contraction index a. The 
constants d^a are fixed according to the contraction pat- 
tern. The basis functions {x;^(r),/i = 1, ...,mb} are also 
chosen to be normalized and Rg is the atomic position in 
the unit cell from where the lobe functions are displaced. 
In the basis set the higher angular momentum terms 
p,d, ... are treated by linear combinations of s- functions 
(lobes) with its origins displaced from the nuclear posi- 
tion by vectors A^q,. (for further details see Ref. p2|). 
For instance, a orbital can be written as a difference 
between two lobe functions Up^ = {ui — U2) so that 

the exponent in the basis set will be 

- R,|2 ^ [r, + (-1)'^ - R,,f ^ (r, - Rs^f 



i=y,z 



(29) 



where A is the magnitude of the displacement of the lobe 
centers, and Rgi is the component of the position vector 



(25) 
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of nucleus s. According to Eq. (18), the Slater potential 
in terms of gaussian basis sets and of the self-consistently 
determined density matrix can be written 



Slater 



1 



2p(r) V 

n.n'n — — oo 



EpO,n 



(30) 



^P<""x°(r)x^'(r) j &^}^-^dv'. 



While the integration with respect to the primed r co- 
ordinate can be carried out analytically [g3[, the second 
integration to calculate the matrix elements according 
to Eq. (|l^) can be accomplished only numerically which 
makes the procedure somewhat time consuming. How- 
ever, further attention will be devoted to simplifying this 
integral using auxiliary fit functions . 

Since we are dealing with a Hartree-Fock-Roothan-Hall 
SCF-LCAO-MO approximation in a basis set represen- 
tation [T^ ], one has to apply the Pock matrix formalism 
and has to form the corresponding Pock matrix elements: 



H 



(31) 



where ii^^,J^^ and FJ][f are the one-electronic. Coulomb 
and exchange contributions to the Pock matrix FJ]^ with 
n unit cell so that 



H 



1 



and 



= {xi\vHir)\x:), 



K'^^ = {xl\vAr)\x:)- 



(32) 



(33) 



(34) 



According to Eq. (|l|) and Eq. (gsj) the stationary one- 
electron periodic wavefunctions {ui{k,r)} are obtained 
in the usual way by solving the eigenvalue equations 

^c,^(fc)j 2 e^'"[F«,-e,(fc)S;:j 1=0, (35) 



with 



s;:. = / dvxiwA^). (36) 



at a finite number of fc points in the Brillouin zone (BZ). 
This gives the band structure ej(fc) and the eigenvectors. 

4. Computational details 

On the basis of the formalism described above the total 
energy and the electronic band structure have been cal- 
culated for polyethylene by various exchange-only meth- 
ods. All of these schemes can be considered as orbital 



dependent methods and results are compared with other 
density functional exchange-only methods like the simple 
exchange-only LDA iXa) procedure. Hartree-Fock cal- 
culations have been carried out as well to make valuable 
comparison with the previous methods. For all the cal- 
culations a new code is used which has been developed in 
our group and the results are compared with other poly- 
mer codes like Mintmire's program jlj]. This new code 
is based on the Erlangen periodic Hartree-Fock program 
[ psf and has been developed by S. Kurth and P. Siile [ p6[ . 
During the calculations the number of fc-points and the 
convergence criteria are set to 25 and 10~^ (the change 
of the density matrix elements from one iteration to the 
next), respectively. For polyethylene the experimental 
" zig-zag" geometry is used JT^ . Since the various proper- 
ties of polyethylene depend significantly on the geometry, 
we give the structure we have used for the calculations: 
the C-C and C-H bond lengths and the CCC, HCH, HCC 
bond angles are set to 1.54A, l.lOA and 113, 108 and 113 
degrees, respectively. Recent geometry optimization cal- 
culations with various DPT functionals provided values 
in very good agreement with these experimental bond- 
ing parameters Three types of basis sets are used: 
dementi's minimal (7S/3P), the larger double-C (9S/5P) 
and (11S/7P) basis sets for carbon and the 4S, 6S and 8S 
for the hydrogen atom, respectively |^^. To make valu- 
able comparison with previously published atomic and 
molecular results j6[0] and also to check the reliability of 
our code we have made calculations with a translation 
vector of about 20.0 a.u. and with the minimal num- 
ber of neighbouring interactions (neig = 1). Using these 
parameters the computed properties must be very close 
to those obtained by atomic or molecular codes. The 
test provided nice agreement with atomic calculations 
obtained by Hartree-Fock, Slater's or SSP exchange-only 
method which confirms the reliability of our code. 

5. Results and discussion 



The comparison is carried out for Be and for LiH [|7|Jlq] 
(Table I) in the above-mentioned atomic and molecular 
limit with our polymer code |^^. The deviation of the 
total energies from HP can be compared with results ob- 
tained by others [Q: we found 12, 5 mHartree for Slater 
and SSP, while Gritsenko et al. obtained 11 and 0.0, re- 
spectively. For LiH the deviations: 0.0 , 0.0 mH, while 
Grabo et al. Q found 6 mH for Slater. The discrepancy 
between our and other's results can be considered rather 
small and is probably due to the different basis set ap- 
plied here. The total energies for C2H4 are compared as 
well in the molecular limit and the following values are 
obtained with Hartree-Fock, Slater and SSP exchange- 
only methods: -77.61066, -77.58735 and -77.60790 a.u. 
One can see that the deviation from HP is 24 mH for 
Slater and only 3 mH for SSP. For the highest occupied 
orbital energies we got -.4273, -.4878 and -.3854 a.u. by 
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HF, Slater and SSP. For LiH the corresponding numbers 
are as follows: -.2914, -.3228 and -.3088 a.u. which are 
compared with Grabo's values 0, -.3017 (HF), -.3150 
(Slater) and -.3011 (KLI). We also compared our cal- 
culated HF total energy and Fermi level with those of 
obtained by other ab initio codes like the Gaussian 94 
package for LiH using the ST0-3G minimal basis 
set. The agreement is acceptable which also indicates 
that our code works in a proper way. 

Our results for some physical properties of polyethy- 
lene are shown in Table H and III. The results are listed 
in terms of the number of the neighbouring interac- 
tions. Actually, it turns out that at least 5 neighbours 
are necessary to get converged results for all the physi- 
cal properties we are interested in. For comparison we 
have also calculated the properties for 8 neighbouring 
cells. In Table HI values are given obtained from Mint- 
mire's polymer code 1 14 as well, which is using the multi- 
ple expansion technique for the neighbouring interaction 
and can therefore be considered as an infinite order ap- 
proach mM. However, in this code a different basis set is 
used (72T/51/1*) which makes the comparison somewhat 
more difficult. 

The comparison of total energies shows that Hartree- 
Fock provides the lowest energies while the various Kohn- 
Sham schemes result in somewhat higher energies. OEP 
leads to total energy values that are upper bounds to 
the HF results and generally the inequality, Ehf < 
EoEP < Ekli < Elsdx holds Q. Actually the differ- 
ence between HF and SSP is 10 mH. For small diatomic 
molecules T. Grabo and E. K. U. Gross Q found a dif- 
ference of about 10 mH between KLI and HF as well, 
which certainly indicates that the SSP potential given 
by Eq. (|2^) is not a bad alternative for KLI. For a larger 
basis set the SSP and HF total energies differ more sig- 
nificantly (14 mH), which is however, much smaller than 
the corresponding difference for Slater (44 mH). We be- 
lieve that the bulk part of the difference is not due to the 
approximation of the exchange potentials applied here, 
but to the different nature of Hartree-Fock and DFT ap- 
proaches in accordance with earlier studies on molecules 
However, in a polymer chain calculation we get 
for the one electronic energies quite significant deviation 
from the Hartree-Fock values (Table II-III). There is a 
trend in these results: the employment of exchange po- 
tentials results in the raising of the Fermi level and the 
falling of the first virtual levels giving a smaller gap when 
compared with Hartree-Fock for polyethylene, re- 
mains almost constant with respect to the increase of the 
number of neighbours taken into account in the neigh- 
bouring approach. 

Neglect of the response part in Eq. ( |l3|) leads to some- 
what higher energies for Etot- Thus the addition of a 
repulsive response part to the attractive bare Slater po- 
tential brings Etot and E^ closer to the corresponding 
Hartree-Fock values ||l^ . The error in Etot obtained with 
the naked Slater potential increases with the basis set 
and reaches 20mH for the largest basis set when com- 



pared with SSP. 

Table II-III represent exchange E^ energies as well. 
Slater yields too negative E^ because of the attractive 
character of v^. The addition of the approximate re- 
pulsive step potential (SSP) brings the E^ values much 
closer to E^^ but still remains somewhat more nega- 
tive (Table III). Gritsenko et al. ||l^ found, however, the 
overcompensation of E^^^ values compared with E^^^ 
for atoms. The inequality E^^^ < E^^ Q is clearly 
reproduced by Slater and SSP with the larger basis set. 

The most striking feature to be mentioned is that the 
highest occupied orbital energies e^n obtained by SSP 
differ from HF values significantly, and are less nega- 
tive than e^"- The values e^"'*'^^ are within the range 
of [e^^;e^"]. These results are rather surprising since 
atomic and molecular calculations show that the val- 
ues obtained by HF or OEP are close to each other 
■ While this comes as a surprise we consider the rais- 
ing of the Fermi level due to the periodic effects, which 
come into play only, however, when one employs local- 
ized exchange-potentials for infinite systems. One can see 
how the periodic effect comes into play when the change 
of Cm is examined in terms of the order of the neighbour- 
ing interactions in Table II. As may be read from Table 
II, when neig = 1 (which case is closer to the molecular 
limit than to a periodic system) , em obtained by Slater 
or SSP is rather close to the HF value, however, when 
neig is increased the Fermi level shifted upward and is 
getting closer to zero progressively. This phenomenon 
is clearly the manifestation of periodic effect on highest 
occupied orbitals. Hartree-Fock does not produce such 
a phenomenon since it exhibits non-local orbital depen- 
dent exchange-potential, which is always deeper than the 
local counterpart, and thus is keeping the Fermi level at 
a deeper energy level. Table II also reports that with one 
neighbour (neig =1) the KS-based methods exhibit too 
wide gap and as the periodicity is building up progres- 
sively, the gap is lowered significantly. 

The Xa results differ more significantly from all the 
other methods, yielding higher total energies and higher 
Fermi levels which is due to the wrong exponential 
long range behaviour of the LDA exchange potential. 
The smallest HOMO-LUMO transition, the fundamen- 
tal (eigenvalue) band gap is found with all the methods 
at the edge of the Brillouin-zone {k = 1). As can be 
seen in Table III, increasing the basis set leads to a de- 
crease of the calculated gap for all the methods, which 
is mainly due to the low first virtual levels when com- 
pared with HF. While the HF gap is still too high the 
calculated Slater's, SSP gaps are too low. provides 
value surprisingly close to the experimental one um while 
in the literature Xa is known for typically giving gaps 
which are too small It is worth to note, however, 
that the HF gap goes through significant changes with 
the size of basis set as well in the [14, 24] eV range 
It is mainly due the downwards movement of the HF first 
virtual level em+i while the Cm level remains almost un- 
changed with respect to the size of the basis set. As is 
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well-known, Hartree-Fock provides too large gaps p3[ . 
Although Slater and SSP give values around the experi- 
mental one (~ 8.8ey p^ , p4 |) with the minimal basis set 
(Table II) the gap obtained by SSP is closer to the Xa 
value. However, the band gap is narrowed significantly 
when the basis set is enlarged (Table III). This clearly 
must be attributed to the downwards movement of the 
first virtual levels tm+i and at the same time to the rais- 
ing of the Fermi level tm- 

To understand the "small eigenvalue gap" problem in 
Fig. 1 and 2 the band structure is plotted obtained 
by various methods using the dementi's minimal and 
double-C bases (Fig. 1 and 2, respectively). Analyzing 
and comparing systematicaly the band structures in Fig 
1-2 one can see that "new" virtual bands appear below 
the first virtual level obtained with the minimal basis 
set. The larger the basis set is, the larger the virtual 
space becomes which can result in the appearance of vir- 
tual states with low orbital energies. We would like to 
emphasize that the appearance of new low-lying virtual 
levels below the ones of the minimal-basis calculations is 
certainly not unexpected. As the basis set increases new 
states appear among the virtual canonical orbitals with 
orbital energies that form a dense subset of possible vir- 
tual energies Q . It is not easy then to separate out the 
virtual state which corresponds to the lowest conduction 
band state and also the one-electron excited states are 
hidden in the virtual subspace. 

A more serious aspect of the small gap problem to 
be considered is that one has to take into account the 
derivative discontinuity to get the right gap when using 
OEP For continuum approximations to E^c hke 

LDA or GGA this discontinuity vanishes. In OEP and 
in our approaches, however, we have a finite derivative 
discontinuity (|,^. In OEP the exact band gap can be 
written as follows felUSdl, 



A — A-'"''"'^ 4- A 

^ ~ '-^nonmt ' '-^xc 



(37) 



(N)- 



.KS 



(N) 



where N is the number of electrons. One has to add 
the xc-contribution (A is identical with the band gap) 
to the difference of the Fermi (ith level) and first virtual 
KS one-electron energies (level v) which represent the 
highest valence band and the energetically lowest conduc- 
tion band states, respectively. A^^-^j^^ is the Kohn-Sham 
eigenvalue band gap. The exchange-only contribution to 
the exact band gap A^; can be given |33] 



A^(z 



-OEP r 



i;-— [p]\u^) 



— {vi\vi) + 2{vi\ii>), 



(38) 



with v^^ being the nonlocal HF exchange operator con- 
structed from the A^/2 occupied KS orbitals. This for- 
mula is coming from the first order DFT perturbation 
theory and the eigenvalue gap A^J^^^j represents 
the zeroth order term in the perturbation expansion [p5[ . 



We calculated A^; using Slater and SSP methods and the 
results can be seen in Table II-III. Values are also ob- 
tained when the numbers of the neighbours are minimal 
{neig = 1). In Eq. (38) the last terms will vanish when 
neig — > oo. Therefore terms (vilvi) and (viliv) will van- 
ish for systems with periodic boundary conditions in the 
limit of an infinite number of unit cells ||35[| . We have 
studied the convergence of A^, and as far as the minimal 
basis set is concerned we found no significant change in 
the magnitude at neig = 5. The discontinuity is guar- 
anteed to be smaller than the true band gap by its defi- 
nition. This is reproduced by our calculations only with 
small basis set. Calculations with the small basis set re- 
sulted in negative value for A^, which corrects the too 
large KS band gaps in the right way (Table II.). With 
the larger basis sets A^, becomes positive and roughly 
represents the twice of the KS band gap. This is in ac- 
cordance with the finding of Stadele et al. that the 
corrected band gaps are close to the Hartree-Fock gap. 
SSP gives the value of 12.7 eV for the true band gap 
which is really close the HF gap (15 eV). As Stadele et 
al. has pointed out, A^^ equals the difference between the 
HF and the OEP eigenvalue gaps, if we assume that the 
difference between the HF and OEP orbitals is negligible 
p5| . Therefore, it is not surprising, that such an approx- 
imation to exchange-only OEP, like the SSP approach, 
provides band gap very similar to the HF one. 

We got the corrected gap 6.9 eV using the Slater's ap- 
proach which is rather close to the Xa value 6.8 eV when 
the basis (9S/5P) is used. Further increase of the ba- 
sis set (11S/9P) results in even larger gap A for all the 
methods considered in this article. Xa LDA provides a 
gap 7.6 eV close to the experimental value 0. This is 
in accordance with the expectations since a continuum 
functional, like the Xa functional, which averages over 
the discontinuity, should yield a gap that is reasonably 
close to the experimental one. The discrepancy between 
the calculated gaps and the experimental one must be 
partly attributed to the correlation contribution to the 
derivative discontinuity (Ac). Although we would not 
like to go into speculations, however, it is quite probable 
that Ac will provide much smaller contribution to the 
true band gap than Aj; . Further attention must be given 
to this in order to determine the magnitude of lS.xc pre- 
cisely. Nevertheless, the naked Slater exchange potential 
provides a corrected band gap A = 9.7 eV in much bet- 
ter agreement with the experiment (~ 8.8 eV) when the 
large (11S/7P) basis set is used (Table HI.). 

The effect of electron correlation on the gaps also 
should be carefully examined, however the recent pub- 
lications indicate only minor importance of correlation 
in DFT gap calculations [^i). This may well be due to 
the small effect of LDA correlation potential on the eigen- 
values H] . Because the correlation energy is strongly de- 
pendent not only on occupied states but on the virtual 
states as well, no simple density functional can be ex- 
pected to yield coherent gaps for different systems. Only 
those correlation energy functionals can provide reason- 
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able eigenvalues and can be checked against experimental 
data by any chance, which are treated strictly at the or- 
bital dependent level of theory. 

The comparison of the Slater and SSP valence energy 
bands with those of Xa and HF again shows that the 
shape of the bands does not change very much (Fig 1- 
2). But the two highest valence energy bands are much 
further apart than in Xa or HF. This is at least in part 
responsible for the small gaps we get for Slater and SSP. 
As already mentioned above, comparing Fig. 1 with Fig. 
2 one can see that using larger basis set (Fig 2.) virtual 
bands are appearing below the zero energy level. This 
phenomenon is clearly demonstrated by all the Kohn- 
Sham based methods except Xa. Note that the lower 
three curves in both figures are the occupied valence en- 
ergy levels (the deepest, core level is not plotted). The 
virtuals of HF do not go below zero as well. However, the 
first virtual levels are lowered significantly as well so that 
the band gap is reduced from 22 eV to 15 eV. The highest 
three virtual levels in Fig. 2 are similar to the virtuals 
obtained by the minimal basis set (Fig. 1). The smallest 
energy difference from the highest occupied level to the 
first virtual is at the edge of the Brillouin zone (fc = 1), 
however, around k — 0.5 one can see quasi-degeneracy 
of certain bands or even crossing of virtual levels close 
to the edge of the Brillouin zone (Fig. 2.). Comparing 
Fig 1 and Fig. 2 it is obvious that increasing the basis 
set shifts the valence bands upwards but does not change 
their shape and relative position very much. But from 
the same figures one can see that the shape and posi- 
tion of the virtual bands changed a lot when increasing 
the basis set. The energy bands get shifted upwards as 
one moves from Slater to SSP. This seems to be consis- 
tent with Gritsenko's and Krieger's work: in Gritsenko's 
Fig. 4 the Slater exchange potential for Ne and Mg 
gets shifted upwards when the response part is added. 
Krieger's Fig. 2 shows the same behaviour for the Ne 
atom as one goes from the Slater exchange potential to 
better exchange potentials as KLI and exact DEP. Com- 
parison of the Slater and SSP valence energy bands with 
those from Xa and HF again show that the shape of the 
bands does not change very much. 

The maximum of the valence band energy is appear- 
ing in the range of [-6,-11] eV with the lowest and highest 
values for Xa and Hartree-Fock, respectively, compared 
to experimentally suggested values for the ionization po- 
tential of 9.6 — 9.8 eV [Q. However, it must be empha- 
sized that others give lower experimental values in the 
7.6 - 8.8 eV range for this property (l|j2|. SSP pro- 
vides the much lower —4.1 eV value (Slater: —7.6 eV). 
By density functional linear muffin-tin orbital method 
the value of —5.1 eV is obtained for helical polyethylene 
p^ . The calculated width of the lowest valence band 
for Xa, HF, SSP and Slater's approach is, respectively 
6.2,9.0, 4.4 and 3.6 eV, which is to be compared with 
the experimental value of 7.2 cV [pl|. The total valence 
bandwidth is 14.0 eV (Xa), 15.3 eV(Slater and SSP) and 
19.8 eV (HF), compared to an experimental value of 16.2 



eV [Q. Slater and SSP perform quite well for this band 
width. The calculated gap between the lowest valence 
band and the minimum of the wider valence bands is 3.3 
eV (HF), 1.9 eV (Xa), 3.0 (SSP) and 3.2 eV (Slater), 
respectively. These numbers have to be compared with 
the experimental 2.0 eV [Q. The HF and Xa values are 
in accordance with those obtained by others |l^ jlj] . The 
bottom of the valence band with a symmetry is roughly 
in the range of [-31,-22] eV, which is in qualitative agree- 
ment with photoemission data [^ with 2 eV difference 
in avarage. 

6. Conclusions 

Calculations have been carried out for polyethylene 
with a new polymer code using various orbital-dependent 
exchange-only potentials. The Slater potential as well as 
its improved version, which incorporates a step-potential 
for the response part of the exact exchange potential, are 
used and the results are compared with those obtained 
by Hartree-Fock or the Xa methods. The addition of the 
step-potential as the response part of the exact exchange 
to the Slater potential results in deeper total energy and 
somewhat wider gap for polyethylene. 

To test the quality and reliability of our code, calcula- 
tions have been performed in the atomic and molecular 
limit. The results agree with those obtained by others. 

In general, we find that the band structure calculated 
by different potentials are similar with little qual- 
itative difference, while all of them differ significantly 
from HF which has valence bands significantly lower, es- 
pecially in the deeper valence region and also exhibits 
higher virtual levels. While the band gap obtained by 
dementi's minimal basis set, is close to the experiment, 
a larger basis leads to gaps which are rather small. How- 
ever, this poor performance for the gap is attached not 
to the approximate nature of the methods applied in this 
paper or not even to OEP but rather to the difficulties 
of gap calculations in extended systems. 

Less negative Fermi levels are calculated by Kohn- 
Sham based orbital dependent methods than by Hartree- 
Fock, which comes as a surprise, because in molecules the 
highest occupied energy levels are very close to those of 
HF. On the basis of progressive improvement of the num- 
bers of neighbours in the neighbouring approach the bulk 
part of this difference must be attributed to periodic ef- 
fects. However, it is unclear yet whether the em values 
move in the right or wrong direction and also the ques- 
tion whether the periodic effect on em has any physical 
relevance remains unanswered. Slater and SSP provide 
eigenvalue band gaps which are too low. In the discus- 
sion section we are speculating on the possible reasons 
for this unexpected phenomenon. The exchange compo- 
nent to the derivative discontinuity is calculated in order 
to correct the Kohn-Sham band gap. These calculations 
indicate that the discontinuity A^, is roughly twice of the 
Kohn-Sham eigenvalue band gap. The corrected band 
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gaps are close to those of obtained by LDA. On the ba- 
sis of results reported here the approximate exchange- 
only potential SSP seems to provide a true band gap for 
polyethylene which is rather close to the Hartree-Fock 
gap. The simple Slater exchange potential provides band 
gap in much better agreement with the experiment. 
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TABLE I. Comparison of HF, Slater and SSP results for Be 
and LiH obtained with the DIOGENES code using dementi's 
minimal basis set with values from the literature. 









Slater''"^ 


Slater""'^'' 


ggpour 


ggpoth, 


Be i?tot/a.u. 
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-14.573 


-14.556 


-14.551 " 


-14.562 


-14.560 


LiH Etot/si.Vi. 


-7.954 


-7.987 ■= 


-7.954 


-7.981 = 


-7.954 


-7.987 




-7.8629 ° 


-7.8620 










€m/a.u. 


-0.2914 


-0.3017 " 


-0.3228 


-0.3150 " 
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-.2870 


-.2857 ° 











'^with ST0-3G basis 

"with ST0-3G basis using Gaussian94 |^ program pack- 
age 
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TABLE II. Calculated properties of polyethylene by vari- 
ous DFT methods using the minimal basis set 
HF, Slater, SSP and Xa denote the Hartree-Fock, ex- 
rchange-only method with Slater's potential (see Eq. 16), 
•"Slater's potential with orbital dependent step potential 
Eq. ( pc| ) and the Xa exchange-only density functional 
cmethod, respectively. All the properties are in a.u. except 
the vertical HOMO-LUMO eigenvalue gap which is given in 
eV and which is computed at the edge of the Brillouin zone 
(fc = 1). Neig is the number of the neighbours taken into 
account in the neighbouring approach. Etot/u, E^ are the 
calculated total and exchange energies per unit cell. Sm and 
em+i axe the highest occupied and the first virtual energy lev- 
els. Ax is the exchange contribution to the total derivative 
discontinuity. The calculations were carried out on the CH2 
unit cell. 



neig 


HF 


Slater 


SSP 


Xq 


1 Etot/u (a.u.) 


-39.03535 


-39.01888 


-39.02768 


-38.41375 


Ex (a.u.) 


-5.82128 


-5.84774 


-5.79080 


-5.15416 


em (a.u.) 


-.6672 


-.4321 


-.4924 


-.5070 


em+i (a.u.) 


.3978 


-.4083 


-.1084 


-.1535 


gap (eV) 


29.0 


15.1 


16.4 


18.0 


A, (eV) 




-5.7 


-4.4 




5 


-38.88012 


-38.86891 


-38.87278 


-38.30351 




-5.75048 


-5.76397 


-5.7379 


-5.16982 




-.4747 


-.3524 


-.2312 


-.2761 




.3444 


-.0433 


.1310 


.1514 




22.3 


8.4 


9.9 


11.6 






-3.5 


-1.9 




8 


-38.88008 


-38.86890 


-38.87275 


-38.30348 




-5.75049 


-5.76400 


-5.7378 


-5.16983 




-.4760 


-.3518 


-.2318 


-.2769 




.3436 


-.0425 


.1303 


.1510 




22.3 


8.4 


9.9 


11.6 






-3.5 


-1.9 





Fig. 1 The calculated valence band structure (eV) 
obtained by various exchange-only methods as a function 
of the dimensionless k variable with fc = being the 
zone center and k = 1 the zone edge. The Brillouin 
zone is that corresponding to the CH2 unit cell. SoUd 
lines correspond to occupied and dashed Hues to virtual 
levels. The dementi's minimal basis set is used and five 
neighbours are considered for the CH2 unit cell. 



Fig. 2 The calculated valence band structure (eV) 
obtained by various exchange-only methods. Solid lines 
correspond to occupied and dashed lines to virtual levels. 
The two lowest dashed curves of virtual levels are of par- 
ticular interest (see text). The dementi's double-C basis 
set is used and 8 neighbours are considered for the CH2 
unit cell. Notations are the same as in Fig. 1 
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TABLE III. Calculated properties of polyethylene by vari- 
ous DFT methods using the dementi's double-zeta basis set 
(9S/5P) and with the (11S/7P) one 



neig 




HF 


Slater 


SSP 


Xa 


5 (9S/5P) 


Etot/n 


-39.01068 


-38.96683 


-38.99704 


-38.41919 




Ex 


-5.8783 


-5.9368 


-5.8841 


-5.2527 




Cm 


-.4049 


-.2951 


-.1628 


-.2227 




em+1 


.1489 


-.1590 


.0007 


.0256 




gap 


15.1 


3.7 


4.4 


6.8 




Ax 




3.4 


3.8 




8 (9S/5P) 




-39.01056 


-38.96684 


-38.99700 


-38.4190 






-5.8780 


-5.9367 


-5.8832 


-5.25221 






-.3984 


-.2774 


-.1490 


-.2116 






.1558 


-.1424 


.0116 


.0375 






15.1 


3.7 


4.4 


6.8 








3.2 


3.7 




8 (11S/7P) 




-39.02326 


-38.99918 


-38.01359 


-38.42701 






-5.8980 


-5.9486 


-5.9036 


-5.2658 






-.3786 


-.2401 


-.1130 


-.2421 






.2320 


-.1151 


.0755 


.0354 






16.6 


3.4 


5.1 


7.6 








6.3 


7.6 




~ oo 


(Mintmire) 








-38.45548 



-5.27507 
-.2111 
.0750 
7.8 



gap (Exp) 8.8 " 

The notations are the same as in Table 2 
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